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Magnetically-driven ferroelectric atomic displacements of the order of lO^'^A have been observed 
in orthorhombic (perovskite like) YMnOa by a single-crystal synchrotron x-ray diffraction. The 
refined polar structure shows the characteristic bond alternation driven by the exchange striction 
in staggered Mn-O-Mn arrays with tt-Ul type ordering, giving rise to a spontaneous polarization 
along a-axis. First-principles calculations based on the Berry phase method as well as on the 
experimentally refined crystal structure can reproduce the observed polarization value. 

PACS numbers: 75.80.-|-q, 77.80.-e, 61.05.cp 



I. INTRODUCTION 

A strong coupling among spin, charge, and lattice de- 
grees of freedom is a key requirement for future spintronic 
devices. Magnetoelectric multiferroics is one such exam- 
ple satisfying this demandii^, as typified by the magnetic- 
order-induced polarization in orthorhombic (perovskite 
like) i?Mn03 (R: Tb and Dy) and NiaVaOg with a cy- 
cloidal spin order'^r^^. The relation between cycloidal spin 
structure and electric polarization has been extensively 
discussed in terms of spin current- or by the inverse 
Dzyaloshinskii-Moriya model-"^, both of which consider 
the spin-orbit coupling as the relevant interaction. Sub- 
sequently, first-principles calculations using the density 
functional theory (DFT) were performed to explain the 
actual value of the spontaneous polarization in a quanti- 
tative manner-iSiii . 

In displacement-type ferroelectrics, where the electric 
polarization had been attributed to the displacements of 
the constituent ions in early days, an essential contri- 
bution of the quantum Berry phase of valence electrons 
was revealed by accurate structural analysis and first- 
principles calculatioE^i2r2i. Therefore, also in the multi- 
ferroic system, the collaboration between the theoretical 
calculation and the crystal structural analysis will give 
critical information to the microscopic origin of the fer- 
roelectric polarization. Nevertheless, little has been clar- 
ified for the crystal structures of the ordered phases in 
multiferroics, primarily due to the extremely small lat- 
tice displacements driven by the magnetic order as re- 
flected in the small polarizations, typically <1 /iC/cm^. 
For a typical multiferroic TbMnOa, the intensities of 
the x-ray superlattice reflections with a wavenumber 2qm 



{qm' wavenumber of spin modulation), which appear de- 
pending on the type of spin helixi^, are four orders of 
magnitude smaller than that of the strong fundamental 
reflections. These are contrastive to the cases of con- 
ventional displacement-type ferroelectrics BaTiOs and 
PbTiOa with large polarization (several tens of /iC/cm^). 

Among the origins of spin-driven ferroelectricity, lat- 
tice striction induced by the symmetric spin exchange 
Si ■ Sj has also been of great interest because of the po- 
tentially large polarizationiii^, which can be ascribed to 
the larger energy scale of the exchange interaction ir- 
relevant to Dzyaloshinskii-Moriya interaction (spin-orbit 
coupling) . Typical examples of exchange-striction-driven 
ferroelectrics are orthorhombic (perovskite like) RMnOs 
{R= Y and Ho-Lu) showing the tt4-4--type (-E-type) spin 
orderly Ca3(Co,Mn)209ia, and GdFeOsJ^. As for the 
orthorhombic i?Mn03, the existence of the enhanced po- 
larization of about 6 /iC/cm^ has been proposed in the 
-E-type phase by the first-principles calculations based 
on the Berry phase method and theoretically optimized 
crystal structur o^°'^^ . Recently, the polarization in the 
-E-type phase has been confirmed for polycrystalline sam- 
ples to be much larger than those of the helimagnetic 
phase, but the value appears one order of magnitude 
smaller than the calculated valu o^^i^^ . Such discrepancy 
should be solved by knowing the accurate crystal struc- 
ture of the multiferroic state as well as by comparing the 
experimental values with the result of the first-principles 
band-theoretical calculation. In this context, the or- 
thorhombic RMnOs with the E-iype order is an ideal sys- 
tem to study the microscopic origin of the magnetically- 
driven ferroelectricity from the structural viewpoint; the 
advantages are that (i) the lattice displacements are ex- 



pected to be relatively large because the polarization is 
particularly large among the magnetically-driven multi- 
ferroics, (ii) extensive reports regarding to the multifer- 
roic properties have been accumulated for the series of the 
orthorhombic KMnO^. In this paper, we report the ob- 
servation of the atomic displacements in the ferroelectric 
phase of orthorhombic YMnOa by means of synchrotron 
x-ray diffraction for the single crystalline sample. The 
atomic displacements with alternating Mn-O-Mn bond 
angles were found to produce polarization along the a- 
axis (in the Phnm setting) in the i?-type phase. By a 
first-principles calculation, we can quantitatively repro- 
duce the experimentally observed polarization value for 
appropriate values of on-site Coulomb interaction U . 



II. EXPERIMENTAL AND COMPUTATIONAL 
PROCEDURES 

Single crystals of orthorhombic YMnOa were synthe- 
sized by a high-pressure treatment (1573 K and 5.5 GPa) 
on the precursor sample of hexagonal YMnOs with us- 
ing a cubic- anvil- type high-pressure apparatus. Details 
of the crystal growth of orthorhombic YMnOs under high 
pressure will be reported elsewhere^^. Although the as- 
prepared crystals with a dimension up to 500 ^m were 
twinned with alternating the a- and 6-axes, the typical 
single-domain size was as large as 100 ^ 200 /im. While 
measurements of the polarization and the dielectric per- 
mittivity were carried out on a twinned crystal, the mag- 
netic susceptibility, infrared reflectivity and synchrotron 
x-ray diffraction were measured for twin-free small crys- 
tals. 

For measurements of the dielectric permittivity and 
the electric polarization, gold electrodes were deposited 
on both faces of the rectangular crystal with dimensions 
of 200 X 240 X 200 iivo? . As a pohng procedure, an elec- 
tric field of up to 12 kV/cm was applied at 40 K, followed 
by cooling to 2 K, and then the electric field was turned 
off. Displacement current was measured with increasing 
temperature at a rate of 5 K/min, and was integrated 
as a function of time to obtain the polarization. For the 
E J- c measurement, electric field was perpendicular to 
the c-axis and slanted at an angle of approximately 33.5° 
from the a- (6-) axis for the major (minor) twin domain; 
this was confirmed by a polarizing microscope. (The ac- 
tual twin domain boundaries are observed along [1 1 0] 
or [1 -1 0] directions.) The reduction of the polarization 
value due to this misorientation is corrected by multi- 
plying the observed raw data by 1.2 (= 1/ cos (27r ggp„ )). 
Magnetic susceptibility was measured in a heating run at 
a magnetic field of 1 kOe with using a SQUID (Super- 
conducting Quantum Interference Device) magnetome- 
ter. A Fourier-transform spectrometer (Bruker IF66v/S) 
combined with a microscope (Hyperion) and a HgCdTe 
detector was used for infrared spectroscopy. 

Synchrotron x-ray diffraction experiments were per- 
formed on BL02B1 at SPring-8, Japan. The photon en- 



ergy of the incident x-rays was tuned at 34.91 keV. X-ray 
beams were shaped into a square of 300x300 /im^ by a 
collimator, which was large enough to irradiate the whole 
sample. A single crystal was crushed into cubes with a 
typical dimension of about 50 /im. The absorption co- 
efficient /i for 34.91 kcV is calculated to be 47.4 cm~^. 
Therefore, ^r was about 0.25, which certified that the 
absorption effect was small enough for the empirical ab- 
sorption corrections^. A large cylindrical imaging plate 
was utilized to measure diffracted intensities. The tem- 
perature was controlled by a helium gas stream cryostat. 
Rapid- Auto program (Rigaku Corp. and MSC.) was used 
to obtain an F-table. Diffraction data were collected up 
to the resolution of 0.2 A. SIR200426 and CRYSTAL- 
STRUCTURE (Rigaku Corp. and MSC.) programs were 
used for analyzing the crystal structure from the F-table. 

As for first-principles calculations, we performed DFT 
simulations using the VASP code^^ within the GGA+U— 
formalism with various U values. In addition, we used the 
Heyd-Scuseria-Ernzerhof (HSE) screened hybrid func- 
tional method, ^^ which mixes the exact non-local Fock 
exchange and the density-functional parametrized ex- 
change. The HSE is known to improve the evaluation 
of the Jahn- Teller distortion, with respect to GGA-|-?7 
approaches?-. The cut-off energy for the plane-wave ex- 
pansion of the wavefunctions was set to 400 eV and a 
fc-point shell of (2, 3, 4) was used for the Brillouin zone 
integration according to Monkhorst-Pack special point 
mesh. Spin-orbit interaction is not taken into account. 



III. RESULTS AND DISCUSSION 

A. Magnetic and dielectric properties 

The temperature dependence of various physical prop- 
erties of YMnOs measured on the single crystal is shown 
in Fig. [TJ Two kink-like anomalies are clearly observed at 
T]<ii—'i2 K and rN2=30 K, respectively, in magnetization 
along the 6-axis, while they are barely discerned along the 
a- and c-axes (Fig. [Ija)), corresponding to the magnetic 
phase transitions from the paramagnetic to the sinusoidal 
and then to the -E-type structures with spins pointing al- 
most along the 6-axisii22ii^. In the E-type phase, the 
ferroelectric polarization perpendicular to the c-axis ap- 
pears, as accompanied with a discontinuous jump of the 
dielectric permittivity e at Tn2, as shown in Figs. HJb) 
and (c). In contrast, the anomaly of e along the c-axis is 
hardly discernible and the polarization along the c-axis is 
very small. These results indicate that ferroelectric po- 
larization lies within the a6-plane in the _E-type phase. 
While the polarization value at 2 K is not saturated at a 
poling field of 12 kV/cm as shown in the inset, the sample 
showed a breakdown for higher poling fields. Neverthe- 
less, the genuine value of polarization at 2 K is expected 
to be at least 0.24 /zC/cm^. 
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FIG. 1: (Color online) The temperature dependence of (a) 
magnetic susceptibility, (b) polarization (P), (c) electric per- 
mittivity (e), (d) spectral weight of the activated 77 meV oxy- 
gen stretching phonon mode (see also Fig. [2](b)), and (e) inte- 
grated intensity of synchrotron x-ray diffraction of a forbidden 
(0 112) reflection and satellites (0 11±S 2). The inset to (b) 
shows the poling electric-field {E) dependence of the P at 2 
K. All the measurements were performed in heating process. 
The temperature dependence of polarization shows a step-like 
behavior which would probably arise from the phase coexis- 
tence occurring near the first-order phase boundary with the 
incommensurate cycloidal phase^^'^'^. The gray lines in (d) 
and (e) are the guide to the eyes. 



B. Infrared reflectivity 

Figure [JJa) shows far-infrared reflectivity spectra of a 
polycrystalline sample. A lot of phonon modes are ob- 
served due to the orthorhombic distortion of the lattice 
from the cubic perovskite structure, where three infrared- 
active phonon modes, i. e. stretching, bending and ex- 
ternal modes, are identified. A clear anomaly is observed 
around 77 meV below rN2 as indicated by an arrow. Ac- 
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FIG. 2: (Color online) Infrared reflectivity spectra of the oxy- 
gen stretching phonon of (a) polycrystalline and single crystal 
samples for (b) _E" || a and (c) E" \\ b both in the -E-type and 
in the sinusoidal phase. Phonon anomaly is observed only 
for E'^ II a configuration. A dip structure around 77 meV 
is attributed to an originally infrared-inactive but lowered- 
lattice-symmetry induced phonon mode. In the inset, the in- 
frared data at 9 K and the fitting spectra with three Lorentz 
oscillators are shown. 



cording to Ref. |32| . a Raman-active mode at 77 meV gains 
the infrared activity due to the lowering of the symmetry 
in the ferroelectric phase. By using an as-grown single- 
crystal surface with a dimension of 50x100 /im^, the po- 
larization of each phonon mode was clearly decomposed, 
as shown in Figs. [2Ib) and (c). The anomaly at 77 meV 
is observed only for S" || a below Tn2 indicating that the 
new optical phonon (oxygen stretching) mode is related 
to the loss of glide-plane normal to the a-axis. There- 
fore, the spectral weight of the mode can be a measure of 
the magnitude of the ferroelectric lattice displacement. 
To deduce the temperature dependence of the spectral 
weight, we fitted the reflectivity spectra above 65 meV 
for the E'^ \\ a configuration by three Lorentz oscillators 



described as S/{uj 



-iujj) (see the inset of Fig. [DJb)). 



Here, S, ujq, and 7 are spectral weight, eigenfrequency, 
and inverse of lifetime of the oscillator, respectively. Be- 
cause of the ambiguity of the absolute value of the reflec- 
tivity in this measurement, the scale factor for absolute 
reflectivity is another fitting parameter. The parameters 
(S, wq, and 7) of two phonon modes were fixed to the 
value obtained by the fitting above rN2 • The parameters 
of the Raman-active phonon mode were obtained for each 
temperature below rN2- The intensity of this mode in- 
creases with decreasing temperature below Tn2 and does 
not vary much below around 20 K, as shown in Fig. [Ifd), 
indicating the saturation of the lattice displacement. 
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FIG. 3: (Color online) Synchrotron x-ray oscillation pho- 
tographs in the (a) paramagnetic (50 K), (b) sinusoidal (30 
K), and (c) S-type phase (21 K). At 21 K(r < Tn2), a (0 11 
2) reflection forbidden in the Pbnm symmetry is clearly ob- 
served, indicating the disappearance of the 6-glide operation. 
At 30 K (T'n2 < T < Tni), (0 11 2) reflection splits along 
the 6*-axis, resulting in two incommensurate (0 11±5 2) su- 
perlattice spots with 5=0.08. Comparison between observed 
(-Fobs) and calculated (-Fiai) structure factor (d) at 21 K in the 
ferroelectric phase and (e) at 50 K in the paraelectric phase. 



C. Synchrotron x-ray diffraction and crystal 
structure analysis 

Synchrotron x-ray oscillation photographs for the re- 
spective phases are displayed in Figs, ^a.), (b), and (c). 
At 50 K above Tni, only fundamental Bragg reflections 
for Pbnm symmetry are observed (Fig. El^a)). In the 
sinusoidal phase, two incommensurate reflections of (0 
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FIG. 4: (Color online) Schematics of the ferroelectric atomic 
displacements projected onto (a) the a6-plane and (b) the 
ac-plane in the ferroelectric B-type phase. Dotted and solid 
arrows indicate the directions and the relative atomic dis- 
placements of Mn and O, respectively. Note that the dis- 
placements of O ions are doubled for the purpose of clarify. A 
unit cell is indicated by a broken line, (c) Changes in the Mn- 
O-Mn bond angles and the distance between the neighboring 
Mn ions in the ab-plane. In the ferroelectric phase (21 K), 
the i?-type spin configuration is represented by large arrows. 
Here, we assumed the ferromagnetic (antiferromagnetic) spin 
arrangement for Mn-O-Mn bond with larger (smaller) bond 
angle. 



11±(5 2) appear, as indicated by arrows in Fig. [Sjb) 
((5=0.08 at 30 K). The intensity and (5-value increases 
and decreases with decreasing temperature, respectively. 
As shown in Fig. [SJc), at 21 K in the i?-type phase, 
the two incommensurate reflections merge into a reflec- 
tion (0 11 2), which is forbidden in the original Pbnm 
symmetry^. As shown in Fig. [Ije) , the integrated inten- 
sity of the additional peaks shows a second-order phase- 
transition like temperature dependence in the sinusoidal 
phase with the incommensurate modulation, while ex- 
hibiting a first-order-like sudden variation in intensity 
and 6 value ((5=0.08 to ^=0) upon the transition from 
the sinusoidal to the £'-type phase. 

By using the data sets of synchrotron x-ray diffrac- 
tion, we performed crystal structure analyses in the E- 
type (i.e. ferroelectric) and the paramagnetic (i.e. para- 
electric) phases, respectively. The comparisons between 
observed and calculated structure factors are shown in 
Figs. E^d) and (e). At 21 K in the ferroelectric phase, 
37305 reflections were observed, and 10500 of them were 
independent. 54 parameters were used for the refine- 
ment. The lattice parameters a, b, and c are 5.24560(10) 
A, 5.82980(10) A, and 7.3295(7) A, respectively. Be- 
cause the forbidden reflections of (0 k I): k=odd with 
Pbnm, setting are observed, the plausible space group is 
P2inm (No. 31), a maximal non- isomorphic subgroup 
of Pbnm. The disappearance of 6-glide reflection allows 



TABLE I: Structure parameters of orthorhombic YMnOa at 21 K in the ferroelectric phase (Space group P2\nm (No. 31)). 
The lattice parameters are a=5. 24560(10) A, 6=5.82980(10) A, and c=7. 3295(7) A, respectively. In the tables, x, y, and z are 
the fractional coordinates. Anisotropic atomic displacement parameters are represented as Un, U22, U33, U12, U13, and f/23. 





site 




X 




y 


z 




Yl 


2a 


0. 


61103(7) 




0.83569(3) 


1/2 




Y2 


2a 


0. 


64624(7) 




0.66579(3) 







Mn 


46 


0.12942(10) 




0.74832(4) 


0.24935(5) 




01 


46 





.4273(3) 




0.57846(18) 


0.3044(2) 




02 


2a 





.0168(3) 




0.7118(2) 


1/2 




03 


2a 





.2409(4) 




0.7868(2) 







04 


46 


0.8291(3) 




0.92158(18) 


0.1971(2) 






t/11 (A^) 


U22 (A^) 




t/33 (A^) 


U12 (A^) 


t/l3 (A^) 


C/23 (A^) 


Yl 


0.00094(3) 


0.00169(2) 




0.00156(4) 


-0.00020(3) 








Y2 


0.00201(4) 


0.00169(3) 




0.00142(4) 


0.00007(3) 








Mn 


0.001510(10) 


0.002200(10) 


1 


1.001450(10) 


0.000270(10) 


0.000040(10) 


0.000170(10) 


01 


0.0031(3) 


0.0033(2) 




0.0039(3) 


0.0009(2) 


-0.0002(2) 


-0.00016(16) 


02 


0.0024(4) 


0.0045(3) 




0.0024(3) 


0.0001(3) 








03 


0.0039(4) 


0.0033(2) 




0.0024(3) 


0.0008(3) 








04 


0.0028(2) 


0.0038(2) 




0.0030(2) 


0.0000(2) 


-0.0004(2) 


-0.00037(15) 



TABLE IL Structure parameters of orthorhombic YMnOs at 50 K in the paraelectric phase (Space group Pbnm (No. 62)). 
The lattice parameters are a=5. 24550(10) A, 6=5.82890(10) A, c=7.3313(7) A. 
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a uniform atomic displacement along the a-axis. In the 
structure analysis, the flack parameter^ is estimated as 
0.503(18), indicating that the volume fractions of two 
kinds of ferroelectric domains are nearly equal. The relia- 
bility factors are i?= 3.16%, i?„=2.74%, GOF(Goodness 
of fit)=1.127. At 50 K in the paraelectric phase, 15409 
reflections were observed, and 5668 of them were inde- 
pendent. 29 parameters were used for the refinement. 
The lattice parameters are as follows: a=5. 24550(10) 
A, 6=5.82890(10) A, c=7.3313(7) A. The space group 
is Phnm (No. 62). The reliability factors are R= 2.41%, 
i?^=2.86%, GOF=1.072. The crystal parameters at 21 K 
in the ferroelectric phase and at 50 K in the paraelectric 
phase are summarized in Tables |I] and |lTl respectively. 

The crystal structure analyses indicated the atomic 
displacements of Mn and O ions within the afo-plane and 
the ac-plane in the ferroelectric phase as shown by ar- 
rows in Figs, mja) and (b), respectively. Here, the po- 
sitions of Y ions are fixed. The displacements of Mn 
ions along the a-axis are observed to be uniform, while 
those along the 6- and c-axes are staggered. The shifts of 



O ions along the a-axis are also uniform but opposite to 
those of Mn ions. The components of the O displacement 
along the h- and c-axes entirely cancel out, similarly to 
those of Mn ions. These uniform atomic displacements 
of Mn and O ions along the a-axis and the associated 
variation of electron wave functions generate the ferro- 
electric polarization parallel to the a-axis. The observed 
atomic displacements along the a-axis are as small as 
0.0041(7) A for Mn ion and 0.0020(16) and 0.0030(16) A 
for two inequivalent O sites, respectively. These values 
are one or two orders of magnitude smaller than those 
of conventional displacement-type ferroelectrics, BaTiOs 
and PbTiOa^^. From the point charge model with nom- 
inal ionic charges (Y; 3-|-, Mn; 3-I-, O; 2-) and these val- 
ues of atomic displacement along the a-axis, ferroelectric 
polarization value is calculated to be approximately 0.5 
/iC/cm^, which is in the same order as the experimen- 
tally observed value. Changes of the Mn-O-Mn bond 
angles and the distance between neighboring Mn ions 
are detailed in Fig. [DJc). At 50 K, there is only one 
kind of Mn-O-Mn bond with an angle of 143.96(5)° in 
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FIG. 5: The absolute value of calculated ferroelectric polariza- 
tion in the i?-type phase by using the Berry phase technique. 
Solid circles represent the results obtained from the GGA+U 
method and the horizontal dashed line corresponds to the re- 
sult from HSE calculations. The simulations were performed 
on the basis of the experimentally refined atomic coordinates 
of the ferroelectric-state crystal (at 21 K). 



basal a6-plane. By contrast, at 21 K, there are two in- 
equivalent bonds of Mn-0(1)-Mn and Mn-0(2)-Mn in ab- 
plane, whose bond angles, 144.09(12)° and 143.78(12)°, 
are slightly larger and smaller than that of Mn-O-Mn in 
the paraelectric phase, respectively. In other words, an 
alternating bond structure emerges along the 6-axis in 
the ferroelectric phase, as shown in Fig. IH^a). 

The mechanism of the polarization generation has been 
proposed to be a symmetric-exchange striction that pro- 
vides the gain of magnetic exchange energy. So as to 
minimize the total nearest-neighbor exchange energy me- 
diated by the local ferromagnetic double-exchange inter- 
action, the bond angles for ferromagnetic tt (antiferro- 
magnetic ti) Mn spin arrangement become large (small) , 
because the transfer integral of eg electron increases (de- 
creases) for the larger (smaller) Mn-O-Mn angle2ii2i. 
The result of the present x-ray structure analysis is in ac- 
cord with the basic idea of the exchange-striction model. 



D. Theoretical calculation 

For a more detailed discussion, we performed first- 
principles simulations of the YMnOa electronic structure. 
In this paper, since the optimization of the geometry - 
in terms of atomic positions - is a rather delicate issue 
(see comment below), we stick to the experimentally re- 
fined crystal structure, as now available from the detec- 
tion of atomic displacements of 10""^ A order. The calcu- 
lated ferroelectric polarization Pa using the Berry phase 
mcthod'^^ for the experimental structure points to —a 
direction, and the absolute value \Pa\ is plotted versus 
the Hubbard parameter U in Fig. [S] The \Pa\ value de- 
creases from 1.8 to 0.5 /xC/cm^ when increasing U from 
to 8 eV. The reduction of \Pa\ is mainly due to the 



reduction of hopping, in turn linked to the reduction of 
the Bg-related double-exchange interaction. The polar- 
ization value is rather strongly affected by the U value, 
which is, unfortunately, experimentally not known for 
YMnOa . It is however experimentally reported that U is 
in the range 5~7.5 eV in the case of a similar compound, 
LaMnOg^i^i, so that the \Pa\ value can be considered to 
be 0.5~0.8 iiC/cm^. In order to overcome the problems 
related to the arbitrary choice of the Hubbard U param- 
eter, we performed hybrid-functional (HSE) calculations 
and obtained an even smaller value of 0.25 /iC/cm^. As 
for the comparison with experiments, we remark that 
theoretical \Pa\ values are larger than the observed values 
~0.08 /xC/cm2 at 21 K and -0.24 /zC/cm^ at 2 K. One 
possible origin for this discrepancy is the insufficiency of 
the poling field in the experiment (see the inset of Fig. [1] 
(b)). The corrected value of the ferroelectric polarization 
observed for the polycrystalline sample is approximately 
0.5 /iC/cni^ near 2 K^'^ in agreement with the theoreti- 
cally predicted order of magnitude. This improvement, 
compared to earlier studie o^"'^^ is due to the use of exper- 
imental structure parameters; the bond angle difference 
is experimentally determined to be 0.3°, while previous 
bare-DFT calculations were based on 3° of bond angle 
difference, obtained from the theoretical optimized crys- 
tal structure. Despite some advances in the comparison 
between theory and experiments, several issues remain 
to be solved, in particular related to the geometry opti- 
mization of the crystal structure. Density functional the- 
ory, while capturing Mn-O-Mn bond-angles and Mn-Mn 
bond-lengths alternating features, indicate d^°i^^ opposite 
direction of Mn and O displacements along the a-axis 
as compared with the experimentally refined structure. 
While the magnitude of the relaxation effects depends 
on the U value, the sign is the same for every U, as 
well as when using hybrid-functionals approaches. These 
results are probably related to the overestimate of p-d 
hybridization for (more itinerant) Cg electrons and/or to 
the underestimate of the interaction between (more local- 
ized) t2g spins. It is possible that magnetostrictive effects 
related to t2g states, expected to produce opposite trends 
for the Mn-O-Mn angles (i.e. larger angles for antiparal- 
lel Mn spins), are not properly described within all DFT, 
DFT+U or HSE approaches. Finally, the assumption of 
a fully collinear spin structure might also play a role in 
the theory-experiments comparison. We therefore hope 
that our work will stimulate further theoretical activities 
aimed at quantitatively reproducing the correct experi- 
mental crystal structure and magnetostrictive effects in 
orthorhombic manganites. 



IV. SUMMARY 

In summary, we have performed measurements of 
the magnetic susceptibility, electric polarization, electric 
permitivity, infrared reflectivity, and synchrotron x-ray 
diffraction to investigate the multiferroic nature of or- 



thorhombic YMnOa by using the single crystal. The re- 
sults indicate that the ferroelectric phase transition with 
the polarization along the a-axis takes place concurrently 
with the E-type magnetic phase transition. In the fer- 
roelectric phase, the phonon anomaly and forbidden re- 
flection are respectively observed in the infrared reflec- 
tivity and synchrotron x-ray diffraction, indicating the 
magnetically-induced lattice distortion; the space group 
changes to P2i nm polar structure. The atomic displace- 
ments of the order of 10~^ A are estimated from the crys- 
tal structure analysis in the ferroelectric phase, indicat- 
ing the formation of the bond alternating structure of 
Mn-O-Mn bond in a6-plane. The feature is qualitatively 
consistent with theoretical predictions of the exchange- 
striction model as a microscopic mechanism for gener- 
ating the electric polarization. Further, by the collabo- 
ration between the crystal structure analysis and first- 
principles calculations, we quantitatively explain the ob- 
served macroscopic polarization value with the input of 
experimental information of the microscopic atomic dis- 
placement. This work may pave a path to understanding 



the microscopic origins of the magnetically-driven fer- 
roelectricity from the structural point of view with the 
support of the first-principles calculations, which would 
facilitate further structural studies of other magnetically- 
driven ferroelectrics with small lattice displacements. 
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